Localization and Mobility Gap in Topological Anderson Insulator 
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It has been proposed that disorder may lead to a new type of topological insulator, called topo- 
logical Anderson insulator (TAI). Here we examine the physical origin of this phenomenon. We 
calculate the topological invariants and density of states of disordered model in a super-cell of 2- 
dimensional HgTe/CdTe quantum well. The topologically non-trivial phase is triggered by a band 
touching as the disorder strength increases. The TAI is protected by a mobility gap, in contrast to 
the band gap in conventional quantum spin Hall systems. The mobility gap in the TAI consists of 
a cluster of non-trivial subgaps separated by almost flat and localized bands. 

PACS numbers: 71.23.-k, 73.21.-b, 73.43.Nq, 



I. INTRODUCTION 

Recent development in topological insulators (TI) has 
greatly enhanced our understanding of topological prop- 
erties in condensed matter [J-Q- Band insulators with 
time reversal symmetry can be classified by a Z2 topolog- 
ical invariant 1/ associated with the occupied bands, 1/ = 
for topologically trivial phase and ly = 1 for non-trivial 
phase[j-i3|- In two dimensions (2D), the TI {v = 1) ex- 
hibits quantum spin Hall effect, whose edge currents are 
robust against weak non-magnetic disorder^, H, [^. This 
dissipation-less transport can only be destroyed by ex- 
tremely strong disorder, which drives the system into 
a traditional Anderson insulator [lol The 2D TI 

has been experimentally realized in HgTe/CdTe quan- 
tum wells, where the thickness of the quantum well can 
be varied to tune the system between TI and normal in- 
sulator [il,[i3. 

Recent numerical simulation reveals an interesting new 
phase, called topological Anderson insulator (TAI) [lo| . 
The TAI is a reentrant TI due to disorder: the disor- 
der drives a 2D topologically trivial insulator into a TI 
phase, then back to trivial insulator at strong disorder. 
This is contrary to the general intuition that disorder al- 
ways tends to localize electronic states. This TAI phase 
has since attracted extensive research interests 14-20j|. 
In the original work, this phase was identified from the 
transport properties showing a two-terminal conductance 
plateau 2e^/ h with extremely small fluctuations [lol| . Fur- 
ther numerical studies confirmed that the plateau con- 
ductance in the TAI is contributed from the dissipation- 
less edge states 1J| , which further suggests the topolog- 
ical origin of this phenomena. Theoretical study based 
on the first Born approximation of the disordered Dirac 
fermions proposed that the TAI originates from a band 
touching and subsequent re-ope ning of a topologically 
nontrivial gap driven by disorder [l5|. The band touch- 
ing has been confirmed in the perturbative and numer- 
ical calculations, but it is not sufficient to explain the 
whole region in the TAI. Very recently, there has been a 
phase diagram for the disordered HgTe/CdTe quantum 



spin Hall well, where the quantum spin Hall phase and 
the TAI are connected [2l|. 

In this paper, we study the topological evolution of the 
TAI, and examine the origin of the TAI from a topolog- 
ical point of view. We calculate the band structure and 
the corresponding Z2 invariants as the disorder strength 
increases. Starting with a topologically trivial insulat- 
ing phase, the bulk gap closes due to the disorder, which 
changes the topological invariants of the occupied bands, 
therefore triggering an insulator-TI transition (band in- 
version). As the disorder strength further increases, a 
bulk gap is re-opened. However, the gap value is too 
small (due to large sample size) and too fiuctuating (due 
to the randomness of disorder) to account for a stable TAI 
phase as observed in transport calculations[l3- We shall 
show clear evidences that the TAI phase corresponds to 
a continuous cluster of nontrivial subgaps, rather than 
a single gap. These subgaps are separated by some ex- 
tremely narrow subbands, and survive through size scal- 
ing and random statistics. In other words, in the TAI 
region, a Fermi level falls into a nontrivial subgap with a 
probability close to 1, regardless of sample size and dis- 
order fluctuations. On the other hand, those extremely 
narrow subbands are strongly localized, therefore they do 
not contribute to electronic transport in the thermody- 
namic limit. This novel phase offers a new realization of 
quantum spin Hall states in solids. 

This paper is organized as follows. In section II, we 
describe the model we use. In section III, the general def- 
inition and calculation methods of topological invariants 
are reviewed. In section IV, the ansatz of defining topo- 
logical invariants for disordered systems is introduced. 
The main results are described in sections V and VI. 



II. THE MODEL 

We first briefly revisit the Bloch's description of elec- 
tronic properties. In real space, the electronic Hamilto- 
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nian in a crystal lattice has the general form 

i af3 (ij) a/3 

(1) 

where i,j are the indices of primary unit cells of the lat- 
tice, and a, /? are the indices of freedom degree within 
the unit cell, e.g., sublattices, orbitals and spins etc. Af- 
ter Fourier transformation Qq, = Sfc '^fca^^'^^S the 
Hamiltonian can be written as 



k ap 



(2) 



where k is defined in the first Brillouin zone. In the 
eigenproblcm 



E 

13 



(3) 



En{k) determines the band structure, and is the 
unit cell periodic part of the Bloch function j'i/'nfc) = 



Un{k)). 



The Bcrncvig-Hughcs-Zhang (BHZ) modell9|, a typi- 
cal tight-binding model with spin-orbit coupling that ex- 
hibits quantum spin Hall phase, is defined on a square 
lattice with one s orbital and one p orbital on each site. 
In the above mentioned representation, the Bloch Hamil- 
tonian is a 4 X 4 matrix written as 



H^p{k) 


f h{k) 


h{k) 


= dol2x2 


9{k) 


f ~l 

^ {a 


do(k) 


= -2D{2- 


rfi(fc) 


= A sin kx , 


dsik) 


= M~2B 



9{k) 
h*{-t 



d2(7y + d^az 



(4) 
(5) 
(6) 



30s kx — COS ky^ 
d2{k) = — Asinfcj, 
I — cos kx — cos ky^ . 

Here a, f3 are the indices of spinorbital within the unit 
cell, a,/? e {1,2,3,4} = {\s t), |p t), k i),\p ;)}. a, 
are Pauli matrices acting on the spinor space spanned by 
s and p orbitals. The real space Hamiltonian Ji of this 
model can be obtained from by a straightforward 
inverse Fourier transformation c^a ~ 7^,- Qae"*'^'^'. 

V y — 

The effect of non-magnetic impurities is included in real 
space by adding a term 



Vi = 



EE 



(7) 



to H, where U{i) are random numbers uniformly dis- 
tributed in (-M^/2,W^/2). 
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1: (Color online) Schematic illustration of Kramers 
pairs. Black is for trivial pair and red for non-trivial pair, 
(a) Band structures E(k). (b) Extension of the Kramers pair 
in (a), represented by the width of a solid bar along the en- 
ergy axis. Different heights (in horizontal direction) of the 
bars are used to distinguish individual pairs. Green bar is for 
nontrivial gap. 



the electronic bands En{k) come in pairs connected at 
time reversal invariant points (TRIPs), which are called 
Kramers pairs (sl [2^. If there are no other degeneracies 
(e.g., disordered "supercells" which will be discussed in 
the following) therefore each Kramers pair (KP) is sep- 
arated from others, a topological invariant can thus be 
defined for each KP[2^. In Fig. [T] (a), we illustrate the 
typical band structures of a time reversal invariant sys- 
tem in the topological aspect. There are 8 bands, forming 
4 KPs, two of which (pairs 2 and 3 in red) are topolog- 
ically nontrivial, which will be defined below. Pairs 3 
and 4 are separated but overlapping in the energy axis. 
We will simply call the gap between them as not fully 
gapped. Most of the information in Fig. [1] (a) can be 
plotted in a simple "bar code" version in [1] (b), where 
the extensions of the KPs and full gaps along the en- 
ergy axis are represented by the width of the bars in this 
direction. 



III. Z2 INVARIANT 



For a time reversal invariant system including both 
spin components, Kramers Theorem states that, all 



In 2D the topological invariant v associated with a KP 
is a Z2 integer defined from the periodic part of the Bloch 
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FIG. 2: (Color online) First Brillouin zone of a time reverse 
symmetric solid in square lattice. Black dots are time reversal 
invariant points. Green region is the effective Brillouin zone 
ri/2, and arrows indicate its boundary dTi/2- Dashed lines: 
mesh used in our calculations (151). 



function u{k) as [22 



— [i dk -A-j dk^F] mod 2, (8) 

27r JdTi/2 -^Ti/a 



A{k)= J2 i{u'{k)\Vk\u'ik)), 
Fik) - fVfc X A(fc) 



(9) 
(10) 



Here rx/2 is the effective Brillouin zone (EBZ) from which 
the rest half can be obtained from its time reverse, dTi/2 
is the boundary of T1/2, as illustrated in Fig. [2] The 
Romantic numbers I and II in equation ([9]) label the 
two branches of a KP (See Fig. [1] (a)). Note the time 
reversal constraint 



\u\-k)) = e\u^\k)) 
|u"(-fe)) = -e\u\k)) 



(11) 



on dTi/2 must be employed for equation (|5]) to make 
sense, where Q = —isy® l2^2K is the time reversal oper- 
ator {sy is the Pauli matrix of physical spin and K here 
is the complex conjugation). A KP is trivial (nontrivial) 
\iv = Q{v = l). The topological invariant of a cluster of 
occupied KPs is just the sum of v of all these KPs in the 
sense of mod 2. Therefore, a gap between two pairs are 
called trivial (nontrivial) if there are even (odd) number 
of nontrivial pairs below it. If a gap is nontrivial, dissi- 
pationless edge states will appear within the ga p, w hen 
the system is truncated with open boundariesjl, 22, 2^ . 
Among several equivalent definitions of Z2 invariant 
[25| , this definition has the advantages of being ex- 
pressed by well known topological quantities, i.e., Berry 



connection A and Berry curvature F 2&[, and being ap- 
propriate for numerical evaluation [23, 27, [2^, which is 
briefly introduced here. After discretizing the EBZ into 
a mesh (dashed lines in Figl2](c)), the field quantities A 
and F can be defined from the eigenstates of the lattice 
sites [2^ [23, [11], based on well-developed lattice gauge 
theories. Note in the numerical calculations, for any mesh 
site k in the EBZ, the phases of the eigenstates (therefore 
the values of the field quantities) are arbitrarily and in- 
dependently determined by the numerical routines {U{1) 
freedom of local gauge choice). Care must be taken to 
cancel all these phase uncertainties when summing up 
the discretized field quantities A and F by equation ([8]) , 
so that the resultant v is gauge independent. Of course, 
the mesh should be dense enough to obtain converged 
values for each KP. 

For the clean systems, the topological properties of this 
model are well understood [Bj, when the lower half bands 
are occupied. When A = 0, the quantum spin Hall phase 
with = 1 is realized when < M/{2B) < 2. When 
tuning M/B, a "band inversion" (ol. [l3| occurs at F point, 
leading to a I-TI transition. The presence of g(fc) breaks 
the conservation of Sz , but the topological invariants does 
not change as long as the finite gap remains. 



IV. ZONE FOLDING 



The topological invariants are defined in k space [22, 
[2^, as introduced above. Impurities break the transla- 
tion invariance of the original lattice and make k badly 
defined. However, for a disordered 2D sample with NxN 
unit cells, the above topological arguments can be re- 
stored if twisted boundary conditions 



TP{r+N-ai) = e*'''"iV(^), ^ir+N-a2) = e 



ikNa.2 



iPir) 
(12) 

are introduced to the opposite boundaries of this fi- 
nite sample, where are primitive vectors of the clean 
lattice [25 



30, 3lj . Physically speaking, this is completely 
equivalent to taking this NxN sample as a large unit 
cell of a 2D super-lattice, so that k can be defined in 
a smaller Brillouin zone with reciprocal vectors bi/N, 
where bi is the reciprocal vector of the original lattice. 
Disorder within this supercell tends to destroy all the 
band degeneracies except those protected by time rever- 
sal symmetry, i.e., Kramers degeneracies. It is reason- 
able to imagine that for sufficiently large N, the topo- 
logical properties of this superlattice can reflect those of 
the "real" disordered system. In the following, we will 
call the primary unit cell of the original clean system as 
a "unit cell" , while the NxN sample as a "supercell" in 
this context. 

In the clean limit, the band structure of the super- 
lattice can be derived directly from that of the original 
lattice by using the standard method of "zone folding". 



4 



which is briefly reviewed here. Now the Bloch Hamilto- 
nian H{k) becomes a AN'^ x AN'^ matrix 



H. 



(13) 



where a, /3 again represent the spinorbital indices within 
the original unit cell, and i,j are the indices of unit cells 
within the supercell. The eigenvalues of H^{k) are re- 
lated with those of the original lattice En{k) as 



I 



N N 
and associated cigenstatcs are 



< < TV - 1 
(14) 



m 
N 



ba) 



V e*(T^''i+T?^''^)-''~x~w„(A; + j^bi + f 62) / 

(15) 

Equations (|14p and (jTSj) can be verified by a straight- 
forward application of Bloch's Theorem, with the new 
definitions of supercell and associated Brillouin zone in 
mind. 



V. A SMALL SUPERCELL 

We will only consider the BHZ model in the case of 
\D\ < \B\, so that the system is always fully gaped be- 
tween the lower and upper halves of bands, when M ^ 0. 
To obtain some insights from analytical treatments, we 
start from a simple stage, a small supercell with 2x2 
unit cells without spin-fiip parts, i.e., A = 0. Now the 
system is decoupled into two sub-systems with single spin 
component, and the topological property of H4x4{k) in 
equation (|4|) can be reduced to that of fe2x2(fc), repre- 

. To further 



sented by spin- resolved Chern number [31 



simplify the analytical treatments, we can only consider 
the spin-up sub-system, because its spin-down counter- 
part can be obtained from a straightforward time reversal 
operation. For this spin-up sub-system with a 2 x 2-site 
supercell, the Hamiltonian is a 8 x 8 matrix 



h^{k) + Vp, 



(16) 



where hfy,g{k) is constructed from the above zone- folding 
technique from the original /i2x2(fe) hi equation ([5|) and 
the impurity term reads 



Vi = W ■ diag(ei, ei, 62, 62, £3, £3, £4, £4) 



(17) 



where are random numbers within the interval 
(1/2,1/2) and W is a single parameter to control the 
disorder strength. This Vj represents random onsite po- 
tential distributed on 4 primary unit cells within the 
2x2 super-cell. We will show that, this minimal model 



in equation (jl6|) that accommodates both disorder and 
topology, can produce some non-trivial results. 

Without impurities, as stated above, the eigenenergies 
and eigcnstates of can be constructed from those of h 
by zone-folding, eqs. and ([15]). The eigenenergies of 
are ordered by their values at the F point fc = as 

-£^1,2..... 8 (r) = -8D- \M -8B\, -AD - |A/-4B|, 

-ad-\m -ab\, -m, m, -ad + \M -AB\, 

-AD+\M-AB\, -8D+\M -8B\, 

with a gap 2M between conductance band E4 and 
valance band E^. The presence of Vi will change band 
structures. Although the band structures including im- 
purities can be deduced from diagonalizing cq. (|16p di- 
rectly, we will treat the impurities as a perturbation, 
therefore the eigenenergies can be expressed as matrix el- 
ements between unperturbed eigenstates \uf). Straight- 
forward calculations show that at F 
order correction to these two states is 



(0,0), the first 



Ei'\T) = (4(F)|F/|4(F)) = ^(ei+e2 + e3 + e4) 

£4), 



W , 



Ei'\r) 



A 

W , 



{ul(T)\Vi\utiT)} = -ie^ 



£2 + £3 



which is just a uniform shift as a simple mean field of 
impurity potentials. The second order correction is 



Ei^'{T) 



where 



E 

1^4 

E 



\{um\Vi\uf{r))[ 

Ei(T)-E,{r) 



-VF^F(ei,e2,e3,£4) 
128{B-D) 



\{uUr)\Vj\ufiT))\' iy2^(ei,62,£3,£4) 



E5(T)-E,ir) 



128{B + D) 



F{ei) = 5(e? + £3 + £3 + £4) - 2(£i£2 + £2£3 



£164 + £364) 



6(£l£3 + £2£4) 



= 0, £1 = £2 = £3 = £4 
> 0, otherwise 



is a scmi-positivc-definite quadrics. Now the gap is 



E„ 



(El + Ei'^ + Ef^) - {El + Ei'> + £;f ^) 



p(2) 



(1) 



2Af 1 



B_ VF^F(el,e2,£3,£4) ^ 
M ' 64(^2 - L>2) ^ 



(18) 



Equation (jl8p is the first important result in this pa- 
per. It suggests that in the weak disorder regime, the 
impurities effectively renormalize m[15| and this renor- 
malization comes from a second order effect of disorder. 
Note the sign of this renormalization term ■ g^gi!fl^i'f 
does not depend on the concrete configuration of random 
impurities. If the clean system is topologically non-trivial 
{M/B > 0), the gap grows as M+const-W'^. This means 
that weak disorder tends to make the two bands expel 
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FIG. 3: Evolution of Kramers pairs and their topological in- 
variants for a 2 X 2 supercell as disorder strength W increase 
(top to bottom), for a definite configuration of {ti}. Bars: 
same as in Fig. [Ub). Model parameters for H in Eqs. (4): 
A = 0.0729, B = -0.0274, C = 0, D = -0.205 = 0.001 
and A = 0, which are also consistent with [lol. [T^. Lattice 
constant is set to be 1. 



each other to avoid a band-touching which wiW trigger a 
transition to a trivial insulator [llj. This is a vivid man- 
ifestation of "robustness against weak disorder" for TI. 
If the clean system is topologically trivial [M / B < 0), 
on the other hand, the gap decays as M — const • W^. If 
M is small, when we tune disorder on, the gap will soon 
close at some small Wc, before strong disorder makes the 
above perturbation treatment unreliable. This gap close 
leads to an I-TI transition with the sign change of M. 
From the topological point of view, the Chern number 
of i?4 changes by 1 after band-touching [nlllilll. Re- 
member we only considered the spin-up block so far, but 
the physical argument of topologically trivial-nontrivial 
transition also applies when the time-reversal invariant H 
involving both spin blocks is considered, by a simple cor- 
respondence between Chern number and Z2 invariant 0, 
as long as these two blocks are decoupled. 

To test the above physical pictures, we calculate the Z2 
topological invariants for a 2 x 2 supercell with both 
spin components included. In Fig. [3l the evolution of 
KPs for a definite configuration of {ei} with increasing 
disorder strength W are plotted. We can see that the 



gap closing predicted by second order perturbation re- 
ally happens at W = 0.024 and it does lead to a topolog- 
ical transition from = to = 1 associated with the 
lower half bands. After this transition, a topologically 
non-trivial gap (the green bar) emerges. This gap will 
develop with further increasing W, until strong disorder 
eventually close it again [llj. This disorder- induced non- 
trivial gap shifts towards positive energy with increasing 
W, as the TAI region observed in [l^l does. This simple 
model itself also paves a route to producing a TI phase 
from a trivial insulator with spin-orbit coupling by con- 
structing a superlattice 



VI. LARGE SUPERCELLS 
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FIG. 4: (a) Two-terminal conductance of a samples with size 
100 X 100, as functions of energy at a given disorder strength 
W = 0.2. Conductance is an average over 300 random con- 
figurations. Conductance plateau with extremely small fluc- 
tuations corresponds to the TAI phase, indicated by the light 
yellow region, (b) The average DOS pave (blue) and the typ- 
ical DOS ptyp (red), calculated from 300 samples with size 
100 X 100 and with periodic boundary conditions. The pa- 
rameters are the same as in Fig. 3 



So far, the origin of TAI seems clear: the disorder trig- 
gers a band touching, or a band inversion, after which 
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a nontrivial gap opens for the TAI phase to live. Un- 
fortunately, this simple argument from a small supercell 
cannot directly be applied to large samples, which will be 
shown in the following. In Fig. 2] (a) , we plot the statis- 
tics of two-terminal conductance of a 100 x 100 samples. 
The two-terminal conductance is calculated by the stan- 
dard non-equilibrium Green's function method (sij. and 
the Fermi energy in the leads are fixed at Ep^'^ = 0.12 to 
offer enough number of channels. The TAI phase is iden- 
tified as the conductance plateau 2e^ /h with extremely 
small fluctuations, as in Ref. [l^. If this region corre- 
sponds to a bulk gap, the density of states (DOS) must 
vanish, at least in the case of periodic boundary condi- 
tion. The single particle local density of states (LDOS) 
is calculated as[37| 



I V =0 |v=1 ~" nontrivial gap Q trivial gap 
-O.OS 0.00 0.05 0.10 0.15 



W=0.02 



The arithmetic mean of the LDOS 



Pavc(£^) =< p{i,E) > 



(20) 



is just the bulk DOS except for a constant factor, where 
^ ■ • • ^ is the arithmetic average over the sites of the 
sample. Meanwhile, the geometric mean of the LDOS 



PtypiE) = exp[< \np{i,E) >] 



(21) 



gives the localization property of the states. In the ther- 
modynamic limit (iV — >■ oo), if Pe,ve(E)/ ptygjE) — > 0, 
then the states around E is localizedjSSl. We thus plot 
Pave and ptyp in Fig- SI (b) in the case of periodic bound- 
ary conditions. Two remarkable features can be read 
from the comparison between Fig. 0] (a) and (b). First, 
the DOS pavG does not vanish in TAI region. As a matter 
of fact, there are regions with smaller DOS outside the 
TAI region. In other words, the TAI phase does not live 
in a bulk gap at all. Second, the TAI region corresponds 
to a vanishing of ptyp. This means that these states are 
extremely localized. 

These surprising results throw a doubt on whether TAI 
can be understood within the above mentioned topolog- 
ical regime. In order to answer this, we repeat the nu- 
merical calculation of v for larger supcrcells. In Fig. [Sj 
we plot the evolution of KPs for a 8 x 8 supercell asso- 
ciated with a definite configuration of {ci}, with increas- 
ing disorder strength W. There is also a band touching 
at W = 0.05, which triggers a nontrivial subgap repre- 
sented by a green bar, as in the case of a small supercell. 
On the other hand, with stronger disorder, for exam- 
ple, around W = 0.2, it develops into a wider region of 
nontrivial subgaps separated by narrow KPs, instead of 
one single nontrivial gap. In Fig. [6j we show the aver- 
age width of KPs within the TAI region, it is clear that 
in the thermodynamic limit, the KPs will be extreniely 
narrow. These narrow KPs are topologically trivial |39l- 
4l| and do not affect the topology (trivial or nontrivial) 
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FIG. 5: Evolution of Kramers pairs and their topological in- 
variants for a 8 X 8 supercell, for a definite configuration of 
{ei}. Disorder strength W increases from top to bottom. Pa- 
rameters are the same as in Fig. 3 



of subgaps between them. In other words, the disorder 
induced nontrivial nature soon hides in the lower KPs 
deeply below the TAI region. We will argue that, those 
narrow and topologically trivial KPs are responsible for 
nonzero DOS in this region, while these nontrivial sub- 
gaps arc responsible for the TAI region observed from 
transport calculations in Ref. [loj . 

The nonzero DOS is easy to understand. Although 
the KPs are microscopically separated, but due to the 
broadening 77 associated with any measure or calculation 
of DOS, they will give rise to a continuous region of finite 
DOS. Moreover, these fiat and well separated KPs tend 
to be strong localized . This is what we have observed 
in Fig. m (b). One may also imagine that a fiuctuating 
transition disorder strength Wc from sample to sample 
might contribute to the statistically nonzero DOS in the 
TAI region. However, as stated in section V, since this 
gap closing is a perturbation effect at weak disorder, the 
ffuctuation of Wc is also very small. Indeed, the numer- 
ical results confirm that (not shown here), compared to 
the width of TAI region, the statistical error of Wc is 
extremely small. 

On the other hand, the origin of TAI phase exhibited 
from transport calculations is more profound. It is well 



7 



0.004 




0) 

> 0.001 
< 



0.000 



FIG. 6: (Color online) Average width of Kramers pairs be- 
tween E — 0.01 and E = 0.03 as a function of supercell size 
A*'. Parameters are the same as in Fig. 3 
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FIG. 7: (Color online) Probability of falling into a nontrivial 
subgap (Q) (lines with symbols) for supercells with different 
size N, at W — 0.2. Every (Q) is averaged over 500 random 
configurations. Green curve is the conductance, same as that 
in Fig. |4] Parameters are the same as in Fig. 3 



known that a nontrivial subgap always gives rise to dis- 
sipationless edge states. However, in our case, these dis- 
order induced nontrivial subgaps are densely and ran- 
domly distributed on the energy axis. To confirm that 
they are indeed responsible for TAX phase, one must ver- 
ify that these nontrivial subgaps can survive through ran- 
dom statistics and size scaling. To characterize this quan- 
titatively, we define a function 

{1, E e a nontrivial and full subgap 
(22) 
0, otherwise 

for a definite supercell size and a definite disorder config- 
uration. The average {Q{E)) over the disorder ensemble 
is the probability for the Fermi energy E to fall into a 



nontrivial subgap. If the above physical picture can make 
sense, {Q{E)) must be a value very close to 1 in the TAI 
region, under random averaging and size scaling. In Fig. 
[7l we show the numerical results of {Q{E)) over averag- 
ing and scaling. The conductance curve is also plotted 
as a comparison. It is happy to see that, the {Q{E)) 
curves converge to a broad peak approaching 1, and that 
the energy region of the broad peak does correspond to 
the conductance plateau of TAI. Note this broad peak 
of {Q{E)) calculated with supercell size > 10 x 10 is 
sufficient to reproduce TAI region identified from con- 
ductance plateau for 100 x 100 samples. In the process 
of scaling [N — > oo), the numbers of subpairs and sub- 
gaps increase, while the widths of individual subpairs and 
subgaps decrease but with different decreasing rate. As 
a result, as disclosed in Fig. [3 for large enough N , the 
total measure of subgaps will dominate over that of flat 
subpairs, {Q) ^ 1 1 — (Q). Fig. [7] is the most impor- 
tant result of this paper. It reveals that, the TAI phase 
corresponds to a cluster of nontrivial subgaps instead of a 
single topologically nontrivial gap. A Fermi energy falls 
into a nontrivial subgap with a probability close to 1. 
The KPs, although contributing to nonzero DOS in this 
region, are so narrow that their measures on the energy 
axis arc extremely small, and they are localized therefore 
do not contribute to the electronic transport. Because of 
the topological origin, it is now clear why the TAI is an 



QSHE phase seen from the real space[21[, and why the 



^ dissipationless currents are still carried by edge states [14 



VII. SUMMARY 

In summary, the topological evolution of TAI is stud- 
ied in a supercell regime. Starting from a trivial insu- 
lator phase with a small gap, weak disorder inevitably 
lead to gap closing between the valence and conduc- 
tion bands, which is a second order perturbation effect. 
This causes an exchange of topological invariants between 
them, and results in a transition to a topologically non- 
trivial phase. In the limit of large supercell, there will 
be very large numbers of subbands and subgaps densely 
distributed on the energy axis. However, there exists 
a continuous region where the Fermi energy falls into 
a nontrivial subgap with an extremely high probability, 
even after a statistical average over the disorder ensem- 
ble. This special region can thus support a stable and 
observable TAI phase. This physical picture also helps 
find disorder-induced topological insulators in other ma- 
terials and higher dimensions (20| . 
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